3.1.48 \(\int \sec ^3(c+d x) (a+a \sin (c+d x))^8 \, dx\) [48]

3.1.48.1 Optimal result
3.1.48.2 Mathematica [A] (verified)
3.1.48.3 Rubi [A] (verified)
3.1.48.4 Maple [A] (verified)
3.1.48.5 Fricas [A] (verification not implemented)
3.1.48.6 Sympy [F(-1)]
3.1.48.7 Maxima [A] (verification not implemented)
3.1.48.8 Giac [B] (verification not implemented)
3.1.48.9 Mupad [B] (verification not implemented)

3.1.48.1 Optimal result

Integrand size = 21, antiderivative size = 121 \[ \int \sec ^3(c+d x) (a+a \sin (c+d x))^8 \, dx=\frac {192 a^8 \log (1-\sin (c+d x))}{d}+\frac {129 a^8 \sin (c+d x)}{d}+\frac {36 a^8 \sin ^2(c+d x)}{d}+\frac {10 a^8 \sin ^3(c+d x)}{d}+\frac {2 a^8 \sin ^4(c+d x)}{d}+\frac {a^8 \sin ^5(c+d x)}{5 d}+\frac {64 a^9}{d (a-a \sin (c+d x))} \]

output
192*a^8*ln(1-sin(d*x+c))/d+129*a^8*sin(d*x+c)/d+36*a^8*sin(d*x+c)^2/d+10*a 
^8*sin(d*x+c)^3/d+2*a^8*sin(d*x+c)^4/d+1/5*a^8*sin(d*x+c)^5/d+64*a^9/d/(a- 
a*sin(d*x+c))
 
3.1.48.2 Mathematica [A] (verified)

Time = 0.16 (sec) , antiderivative size = 111, normalized size of antiderivative = 0.92 \[ \int \sec ^3(c+d x) (a+a \sin (c+d x))^8 \, dx=\frac {a^8 \sec ^2(c+d x) (1-\sin (c+d x)) (1+\sin (c+d x)) \left (192 \log (1-\sin (c+d x))+\frac {64}{1-\sin (c+d x)}+129 \sin (c+d x)+36 \sin ^2(c+d x)+10 \sin ^3(c+d x)+2 \sin ^4(c+d x)+\frac {1}{5} \sin ^5(c+d x)\right )}{d} \]

input
Integrate[Sec[c + d*x]^3*(a + a*Sin[c + d*x])^8,x]
 
output
(a^8*Sec[c + d*x]^2*(1 - Sin[c + d*x])*(1 + Sin[c + d*x])*(192*Log[1 - Sin 
[c + d*x]] + 64/(1 - Sin[c + d*x]) + 129*Sin[c + d*x] + 36*Sin[c + d*x]^2 
+ 10*Sin[c + d*x]^3 + 2*Sin[c + d*x]^4 + Sin[c + d*x]^5/5))/d
 
3.1.48.3 Rubi [A] (verified)

Time = 0.29 (sec) , antiderivative size = 108, normalized size of antiderivative = 0.89, number of steps used = 5, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.190, Rules used = {3042, 3146, 49, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \sec ^3(c+d x) (a \sin (c+d x)+a)^8 \, dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \int \frac {(a \sin (c+d x)+a)^8}{\cos (c+d x)^3}dx\)

\(\Big \downarrow \) 3146

\(\displaystyle \frac {a^3 \int \frac {(\sin (c+d x) a+a)^6}{(a-a \sin (c+d x))^2}d(a \sin (c+d x))}{d}\)

\(\Big \downarrow \) 49

\(\displaystyle \frac {a^3 \int \left (\frac {64 a^6}{(a-a \sin (c+d x))^2}-\frac {192 a^5}{a-a \sin (c+d x)}+\sin ^4(c+d x) a^4+8 \sin ^3(c+d x) a^4+30 \sin ^2(c+d x) a^4+72 \sin (c+d x) a^4+129 a^4\right )d(a \sin (c+d x))}{d}\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {a^3 \left (\frac {64 a^6}{a-a \sin (c+d x)}+\frac {1}{5} a^5 \sin ^5(c+d x)+2 a^5 \sin ^4(c+d x)+10 a^5 \sin ^3(c+d x)+36 a^5 \sin ^2(c+d x)+129 a^5 \sin (c+d x)+192 a^5 \log (a-a \sin (c+d x))\right )}{d}\)

input
Int[Sec[c + d*x]^3*(a + a*Sin[c + d*x])^8,x]
 
output
(a^3*(192*a^5*Log[a - a*Sin[c + d*x]] + 129*a^5*Sin[c + d*x] + 36*a^5*Sin[ 
c + d*x]^2 + 10*a^5*Sin[c + d*x]^3 + 2*a^5*Sin[c + d*x]^4 + (a^5*Sin[c + d 
*x]^5)/5 + (64*a^6)/(a - a*Sin[c + d*x])))/d
 

3.1.48.3.1 Defintions of rubi rules used

rule 49
Int[((a_.) + (b_.)*(x_))^(m_.)*((c_.) + (d_.)*(x_))^(n_.), x_Symbol] :> Int 
[ExpandIntegrand[(a + b*x)^m*(c + d*x)^n, x], x] /; FreeQ[{a, b, c, d}, x] 
&& IGtQ[m, 0] && IGtQ[m + n + 2, 0]
 

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 3042
Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 3146
Int[cos[(e_.) + (f_.)*(x_)]^(p_.)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m 
_.), x_Symbol] :> Simp[1/(b^p*f)   Subst[Int[(a + x)^(m + (p - 1)/2)*(a - x 
)^((p - 1)/2), x], x, b*Sin[e + f*x]], x] /; FreeQ[{a, b, e, f, m}, x] && I 
ntegerQ[(p - 1)/2] && EqQ[a^2 - b^2, 0] && (GeQ[p, -1] ||  !IntegerQ[m + 1/ 
2])
 
3.1.48.4 Maple [A] (verified)

Time = 1.08 (sec) , antiderivative size = 128, normalized size of antiderivative = 1.06

method result size
parallelrisch \(\frac {a^{8} \left (-30720 \left (\sin \left (d x +c \right )-1\right ) \ln \left (\sec ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )+61440 \left (\sin \left (d x +c \right )-1\right ) \ln \left (\tan \left (\frac {d x}{2}+\frac {c}{2}\right )-1\right )+18 \sin \left (5 d x +5 c \right )-8095 \cos \left (2 d x +2 c \right )+166 \cos \left (4 d x +4 c \right )-\cos \left (6 d x +6 c \right )-27580 \sin \left (d x +c \right )-1130 \sin \left (3 d x +3 c \right )+7930\right )}{160 d \left (\sin \left (d x +c \right )-1\right )}\) \(128\)
risch \(-192 i a^{8} x -\frac {1093 i a^{8} {\mathrm e}^{i \left (d x +c \right )}}{16 d}+\frac {1093 i a^{8} {\mathrm e}^{-i \left (d x +c \right )}}{16 d}-\frac {384 i a^{8} c}{d}-\frac {128 i a^{8} {\mathrm e}^{i \left (d x +c \right )}}{\left (-i+{\mathrm e}^{i \left (d x +c \right )}\right )^{2} d}+\frac {384 a^{8} \ln \left (-i+{\mathrm e}^{i \left (d x +c \right )}\right )}{d}+\frac {a^{8} \sin \left (5 d x +5 c \right )}{80 d}+\frac {a^{8} \cos \left (4 d x +4 c \right )}{4 d}-\frac {41 a^{8} \sin \left (3 d x +3 c \right )}{16 d}-\frac {19 a^{8} \cos \left (2 d x +2 c \right )}{d}\) \(176\)
derivativedivides \(\frac {a^{8} \left (\frac {\sin ^{9}\left (d x +c \right )}{2 \cos \left (d x +c \right )^{2}}+\frac {\left (\sin ^{7}\left (d x +c \right )\right )}{2}+\frac {7 \left (\sin ^{5}\left (d x +c \right )\right )}{10}+\frac {7 \left (\sin ^{3}\left (d x +c \right )\right )}{6}+\frac {7 \sin \left (d x +c \right )}{2}-\frac {7 \ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{2}\right )+8 a^{8} \left (\frac {\sin ^{8}\left (d x +c \right )}{2 \cos \left (d x +c \right )^{2}}+\frac {\left (\sin ^{6}\left (d x +c \right )\right )}{2}+\frac {3 \left (\sin ^{4}\left (d x +c \right )\right )}{4}+\frac {3 \left (\sin ^{2}\left (d x +c \right )\right )}{2}+3 \ln \left (\cos \left (d x +c \right )\right )\right )+28 a^{8} \left (\frac {\sin ^{7}\left (d x +c \right )}{2 \cos \left (d x +c \right )^{2}}+\frac {\left (\sin ^{5}\left (d x +c \right )\right )}{2}+\frac {5 \left (\sin ^{3}\left (d x +c \right )\right )}{6}+\frac {5 \sin \left (d x +c \right )}{2}-\frac {5 \ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{2}\right )+56 a^{8} \left (\frac {\sin ^{6}\left (d x +c \right )}{2 \cos \left (d x +c \right )^{2}}+\frac {\left (\sin ^{4}\left (d x +c \right )\right )}{2}+\sin ^{2}\left (d x +c \right )+2 \ln \left (\cos \left (d x +c \right )\right )\right )+70 a^{8} \left (\frac {\sin ^{5}\left (d x +c \right )}{2 \cos \left (d x +c \right )^{2}}+\frac {\left (\sin ^{3}\left (d x +c \right )\right )}{2}+\frac {3 \sin \left (d x +c \right )}{2}-\frac {3 \ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{2}\right )+56 a^{8} \left (\frac {\left (\tan ^{2}\left (d x +c \right )\right )}{2}+\ln \left (\cos \left (d x +c \right )\right )\right )+28 a^{8} \left (\frac {\sin ^{3}\left (d x +c \right )}{2 \cos \left (d x +c \right )^{2}}+\frac {\sin \left (d x +c \right )}{2}-\frac {\ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{2}\right )+\frac {4 a^{8}}{\cos \left (d x +c \right )^{2}}+a^{8} \left (\frac {\sec \left (d x +c \right ) \tan \left (d x +c \right )}{2}+\frac {\ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{2}\right )}{d}\) \(442\)
default \(\frac {a^{8} \left (\frac {\sin ^{9}\left (d x +c \right )}{2 \cos \left (d x +c \right )^{2}}+\frac {\left (\sin ^{7}\left (d x +c \right )\right )}{2}+\frac {7 \left (\sin ^{5}\left (d x +c \right )\right )}{10}+\frac {7 \left (\sin ^{3}\left (d x +c \right )\right )}{6}+\frac {7 \sin \left (d x +c \right )}{2}-\frac {7 \ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{2}\right )+8 a^{8} \left (\frac {\sin ^{8}\left (d x +c \right )}{2 \cos \left (d x +c \right )^{2}}+\frac {\left (\sin ^{6}\left (d x +c \right )\right )}{2}+\frac {3 \left (\sin ^{4}\left (d x +c \right )\right )}{4}+\frac {3 \left (\sin ^{2}\left (d x +c \right )\right )}{2}+3 \ln \left (\cos \left (d x +c \right )\right )\right )+28 a^{8} \left (\frac {\sin ^{7}\left (d x +c \right )}{2 \cos \left (d x +c \right )^{2}}+\frac {\left (\sin ^{5}\left (d x +c \right )\right )}{2}+\frac {5 \left (\sin ^{3}\left (d x +c \right )\right )}{6}+\frac {5 \sin \left (d x +c \right )}{2}-\frac {5 \ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{2}\right )+56 a^{8} \left (\frac {\sin ^{6}\left (d x +c \right )}{2 \cos \left (d x +c \right )^{2}}+\frac {\left (\sin ^{4}\left (d x +c \right )\right )}{2}+\sin ^{2}\left (d x +c \right )+2 \ln \left (\cos \left (d x +c \right )\right )\right )+70 a^{8} \left (\frac {\sin ^{5}\left (d x +c \right )}{2 \cos \left (d x +c \right )^{2}}+\frac {\left (\sin ^{3}\left (d x +c \right )\right )}{2}+\frac {3 \sin \left (d x +c \right )}{2}-\frac {3 \ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{2}\right )+56 a^{8} \left (\frac {\left (\tan ^{2}\left (d x +c \right )\right )}{2}+\ln \left (\cos \left (d x +c \right )\right )\right )+28 a^{8} \left (\frac {\sin ^{3}\left (d x +c \right )}{2 \cos \left (d x +c \right )^{2}}+\frac {\sin \left (d x +c \right )}{2}-\frac {\ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{2}\right )+\frac {4 a^{8}}{\cos \left (d x +c \right )^{2}}+a^{8} \left (\frac {\sec \left (d x +c \right ) \tan \left (d x +c \right )}{2}+\frac {\ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{2}\right )}{d}\) \(442\)

input
int(sec(d*x+c)^3*(a+a*sin(d*x+c))^8,x,method=_RETURNVERBOSE)
 
output
1/160*a^8*(-30720*(sin(d*x+c)-1)*ln(sec(1/2*d*x+1/2*c)^2)+61440*(sin(d*x+c 
)-1)*ln(tan(1/2*d*x+1/2*c)-1)+18*sin(5*d*x+5*c)-8095*cos(2*d*x+2*c)+166*co 
s(4*d*x+4*c)-cos(6*d*x+6*c)-27580*sin(d*x+c)-1130*sin(3*d*x+3*c)+7930)/d/( 
sin(d*x+c)-1)
 
3.1.48.5 Fricas [A] (verification not implemented)

Time = 0.29 (sec) , antiderivative size = 130, normalized size of antiderivative = 1.07 \[ \int \sec ^3(c+d x) (a+a \sin (c+d x))^8 \, dx=-\frac {4 \, a^{8} \cos \left (d x + c\right )^{6} - 172 \, a^{8} \cos \left (d x + c\right )^{4} + 2192 \, a^{8} \cos \left (d x + c\right )^{2} - 1119 \, a^{8} - 3840 \, {\left (a^{8} \sin \left (d x + c\right ) - a^{8}\right )} \log \left (-\sin \left (d x + c\right ) + 1\right ) - {\left (36 \, a^{8} \cos \left (d x + c\right )^{4} - 592 \, a^{8} \cos \left (d x + c\right )^{2} - 2399 \, a^{8}\right )} \sin \left (d x + c\right )}{20 \, {\left (d \sin \left (d x + c\right ) - d\right )}} \]

input
integrate(sec(d*x+c)^3*(a+a*sin(d*x+c))^8,x, algorithm="fricas")
 
output
-1/20*(4*a^8*cos(d*x + c)^6 - 172*a^8*cos(d*x + c)^4 + 2192*a^8*cos(d*x + 
c)^2 - 1119*a^8 - 3840*(a^8*sin(d*x + c) - a^8)*log(-sin(d*x + c) + 1) - ( 
36*a^8*cos(d*x + c)^4 - 592*a^8*cos(d*x + c)^2 - 2399*a^8)*sin(d*x + c))/( 
d*sin(d*x + c) - d)
 
3.1.48.6 Sympy [F(-1)]

Timed out. \[ \int \sec ^3(c+d x) (a+a \sin (c+d x))^8 \, dx=\text {Timed out} \]

input
integrate(sec(d*x+c)**3*(a+a*sin(d*x+c))**8,x)
 
output
Timed out
 
3.1.48.7 Maxima [A] (verification not implemented)

Time = 0.19 (sec) , antiderivative size = 97, normalized size of antiderivative = 0.80 \[ \int \sec ^3(c+d x) (a+a \sin (c+d x))^8 \, dx=\frac {a^{8} \sin \left (d x + c\right )^{5} + 10 \, a^{8} \sin \left (d x + c\right )^{4} + 50 \, a^{8} \sin \left (d x + c\right )^{3} + 180 \, a^{8} \sin \left (d x + c\right )^{2} + 960 \, a^{8} \log \left (\sin \left (d x + c\right ) - 1\right ) + 645 \, a^{8} \sin \left (d x + c\right ) - \frac {320 \, a^{8}}{\sin \left (d x + c\right ) - 1}}{5 \, d} \]

input
integrate(sec(d*x+c)^3*(a+a*sin(d*x+c))^8,x, algorithm="maxima")
 
output
1/5*(a^8*sin(d*x + c)^5 + 10*a^8*sin(d*x + c)^4 + 50*a^8*sin(d*x + c)^3 + 
180*a^8*sin(d*x + c)^2 + 960*a^8*log(sin(d*x + c) - 1) + 645*a^8*sin(d*x + 
 c) - 320*a^8/(sin(d*x + c) - 1))/d
 
3.1.48.8 Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 275 vs. \(2 (120) = 240\).

Time = 0.38 (sec) , antiderivative size = 275, normalized size of antiderivative = 2.27 \[ \int \sec ^3(c+d x) (a+a \sin (c+d x))^8 \, dx=-\frac {2 \, {\left (480 \, a^{8} \log \left (\tan \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )^{2} + 1\right ) - 960 \, a^{8} \log \left ({\left | \tan \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right ) - 1 \right |}\right ) + \frac {160 \, {\left (9 \, a^{8} \tan \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )^{2} - 20 \, a^{8} \tan \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right ) + 9 \, a^{8}\right )}}{{\left (\tan \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right ) - 1\right )}^{2}} - \frac {1096 \, a^{8} \tan \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )^{10} + 645 \, a^{8} \tan \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )^{9} + 5840 \, a^{8} \tan \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )^{8} + 2780 \, a^{8} \tan \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )^{7} + 12120 \, a^{8} \tan \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )^{6} + 4286 \, a^{8} \tan \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )^{5} + 12120 \, a^{8} \tan \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )^{4} + 2780 \, a^{8} \tan \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )^{3} + 5840 \, a^{8} \tan \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )^{2} + 645 \, a^{8} \tan \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right ) + 1096 \, a^{8}}{{\left (\tan \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )^{2} + 1\right )}^{5}}\right )}}{5 \, d} \]

input
integrate(sec(d*x+c)^3*(a+a*sin(d*x+c))^8,x, algorithm="giac")
 
output
-2/5*(480*a^8*log(tan(1/2*d*x + 1/2*c)^2 + 1) - 960*a^8*log(abs(tan(1/2*d* 
x + 1/2*c) - 1)) + 160*(9*a^8*tan(1/2*d*x + 1/2*c)^2 - 20*a^8*tan(1/2*d*x 
+ 1/2*c) + 9*a^8)/(tan(1/2*d*x + 1/2*c) - 1)^2 - (1096*a^8*tan(1/2*d*x + 1 
/2*c)^10 + 645*a^8*tan(1/2*d*x + 1/2*c)^9 + 5840*a^8*tan(1/2*d*x + 1/2*c)^ 
8 + 2780*a^8*tan(1/2*d*x + 1/2*c)^7 + 12120*a^8*tan(1/2*d*x + 1/2*c)^6 + 4 
286*a^8*tan(1/2*d*x + 1/2*c)^5 + 12120*a^8*tan(1/2*d*x + 1/2*c)^4 + 2780*a 
^8*tan(1/2*d*x + 1/2*c)^3 + 5840*a^8*tan(1/2*d*x + 1/2*c)^2 + 645*a^8*tan( 
1/2*d*x + 1/2*c) + 1096*a^8)/(tan(1/2*d*x + 1/2*c)^2 + 1)^5)/d
 
3.1.48.9 Mupad [B] (verification not implemented)

Time = 6.02 (sec) , antiderivative size = 97, normalized size of antiderivative = 0.80 \[ \int \sec ^3(c+d x) (a+a \sin (c+d x))^8 \, dx=\frac {192\,a^8\,\ln \left (\sin \left (c+d\,x\right )-1\right )-\frac {64\,a^8}{\sin \left (c+d\,x\right )-1}+129\,a^8\,\sin \left (c+d\,x\right )+36\,a^8\,{\sin \left (c+d\,x\right )}^2+10\,a^8\,{\sin \left (c+d\,x\right )}^3+2\,a^8\,{\sin \left (c+d\,x\right )}^4+\frac {a^8\,{\sin \left (c+d\,x\right )}^5}{5}}{d} \]

input
int((a + a*sin(c + d*x))^8/cos(c + d*x)^3,x)
 
output
(192*a^8*log(sin(c + d*x) - 1) - (64*a^8)/(sin(c + d*x) - 1) + 129*a^8*sin 
(c + d*x) + 36*a^8*sin(c + d*x)^2 + 10*a^8*sin(c + d*x)^3 + 2*a^8*sin(c + 
d*x)^4 + (a^8*sin(c + d*x)^5)/5)/d